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Abstract. We present linear and circular polarization observations of the ILO masers in 4 distinct regions spread over 1x2 
arcseconds around the HW2 high-mass young stellar object in the Cepheus A star-forming region. We find magnetic fields 
between 100-500 mG in the central maser region, which has been argued to trace a circumstellar disk. The masers further from 
HW2 have field strengths between 30-100 mG. In all cases the magnetic field pressure is found to be similar to the dynamic 
pressure, indicating that the magnetic field is capable of controlling the outflow dynamics around HW2. In addition to several 
H2O maser complexes observed before, we also detect a new maser filament, 1" 690 AU) East of HW2, which we interpret 
as a shocked region between the HW2 outflow and the surrounding medium. We detect a linear polarization gradient along the 
filament as well as a reversal of the magnetic field direction. This is thought to mark the transition between the magnetic field 
associated with the outflow and that found in the surrounding molecular cloud. In addition to the magnetic field we determine 
several other physical properties of the maser region, including density and temperatures as well as the maser beaming angles. 
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1. Introduction 

While the process of low-mass star-formation has been well 
studied, high-mass star-formation is still poorly understood. 
Although several theories propose the formation of high- 
mass stars from the merg er of sever al low-mass young stel- 
' lar objects (e.g. iBonnell et all Il998l) recent studies and ob- 
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servations suggest that high-mass stars form, similar to low- 
mass stars, through accretion from a circumstellar disk ( e.g. 
iMcKee & Tarl l2003l IPatel et all I20051 Ijiang et all Eool . In 
the prevailing picture of low-mass star-formation out of dense 
molecular clouds, strong magnetic fields support the clouds 
against a gravitational collapse. When self-gravity overcomes 
the magnetic pressu re in the cloud core, the formation of pro- 
tostar s ensues (e.g. IShu et al 1 Il987t Mouschovias & CiolekL 
Additionally, magnetic fields likely play an impor- 
tant role in many other stages of star-formation, such as the 
formation of bi-polar outflo ws and a circumstellar disk (e.g. 
lAkeson & Carlstroml Il997l) . Thus, accurate measurements of 
the magnetic field strength and structure in the densest areas of 
star-forming regions (SFRs) are needed to investigate the exact 
role of the magnet ic field in both high- and low-mass starfor- 
mation (see, e.g ISarma et alll200lll2002l) . 

Send offprint requests to: WV (wouter@jb.man.ac.uk) 



Through polarization observations, masers are excellent 
probes of magnetic field strength and structure in masing re- 
gions. For example, polarimetric SiO, H2O and OH maser 
observations in the envelopes of evolved stars have revealed 
the strength and structure of the m agnetic fields during 
the end-stages of stellar evolut i on (e. g Kemball & Diamond, 
Il997t lEtoka & Diamond! |2004 Iviemmings et all l2005t) and 
H2O maser polarization observations have provided stringent 
upper limit s of the magnetic f ield in the megamaser galaxy 
NGC 4258 JModiaz et alll2005h . SFRs also show a rich variety 
of maser species, including OH and H2O . The OH masers are 
often found at several hundred to thousands AU from the SFR 
cores where the density «h, is less than a few times 10 8 crrT 3 . 
Observations of the Zeeman effect on OH masers have been 
use d to determine the SF R magn etic field in those regions (e.g 
ICohen et all Il99fi iBartkiewicz et all l2005h . The H 2 maser 
emission in SFRs is often associated with shocks created by the 
outflows of young stellar obje cts (YSOs) or with a circumstel- 
lar di sk jTorrelles et allll996[ hereafter T96: lGallimore et all 
2003, hereafter G03). The H2O masers are excited in the 
dense parts of SFRs, with number densities «h, be tween ap- 
proximately 10 8 and 10 10 cm -3 felitzur et all Il9 89i). Because 
they are typically small (~1 AU), have a narrow velocity width 
(~ 1 km s _1 ) and have a high brightness temperature T\, > 10 9 K 
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Fig. 1. Total intensity (I) map with polarization vectors of our 
polarization calibrator J2202+4216 (BL Lac). The position an- 
gle (x) of the vectors has been rotated by 77° so that it corre- 
sponds to the VLBA calibration observation. 

(e.g. iReid & Moranl Il98ll) . H2O masers can be used to ex- 
amine the small scale magnetic field strength and structure 
in dense parts of SFRs with polarimetric very long baseline 
interferometry (VLBI) observations. Previous VLBI observa- 
tions have studied the linear polarization of H2O masers as 
tracer of the magnetic field morphology in the SFR s W51 M 
jLeppanen etallT T998). Orion KL and W3 IRS 5 Jlmai et all 
I2003I) . The circular polarization due to Zeeman splitting of the 
22 GH z H2O masers was first observed by Fiebig & Gtistenl 
(1989) with the Effelsberg 100m telescope. These observations 
were confirmed with VLBI bv ISarma et al.l yQOl), who ob- 
served the H2O maser circular polarization in W3 IRS 5 with 
the Ver y Long Basel i ne Ar ray (VLBA). At lower spatial res- 
olution, ISarma et alJ (|2002), also used the Very Large Array 
(VLA) to determine magnetic field strengths in a number of 
SFRs from H2O maser observations. Here we present VLBA 
linear and circular polarization observations of the H2O maser 
structures in the SFR Cepheus A HW2. 

Cephe us A is a high -mass SFR located at a distance of 
~725 pc J.Tohnsorl Il957h . which cont ains a large number of 
radio continuum sources (HW sources: [kughes & Wouterlool 
Il984l) . Additionally it exhibits multi-polar outflows, NH 3 
clouds, Herbig-Haro (HH) objects and infrared sources and a 
complex structure of OH, H2O and methanol masers. The HW 
sources are compact HII regions that are thought to be excited 
by a YSO either externally or embedded in the HII cloud itself 



JCohen et allll984tlGarav et all 1 19961) The brightest of these 
sources is HW2 jRo^rigueTeT^nTl994l) . which is thought to 
contain the main exciting source in the SFR. Surrounding it is 
a rich structure of H? Q masers which has been s tudied in great 
detail (e.g T96: G03: lTorrelles et all fl998l 1200 1 albl hereafter 
T98, TOla and TOlb). More H2O maser structures are found in 
clusters around other H W sources (HW3b and HW3 d), ^ 4-5" 
south from HW2 (T98: ITada et all fl98ll ICohen et all Il984 
iRowland fc Coherl M 984 . The" main, large scale, H2O maser 
structure in the direction of HW2 was interpreted as tracing a 
300 AU radius circumstellar disk perpendicular to the HW2 ra- 
dio jet (T96). Recently, a flattened disk-like structure of dust 
and molecular gas with radius 330 AU oriented perpendicu- 
lar to and spatially coincident with the HW2 ra dio jet has been 
reported (Pate l et allEooilCuriel et alll2005l) . 

Here we examine the polarization properties of the 
H2O masers around Cepheus A HW2 and determine the mag- 
netic field strength and structure. Additionally we describe the 
physical properties of the H2O maser regions and discuss the 
detection of a new H2O maser filament approximately 1" East 
of the HW2 region. 

The observations are described in § [21 and the results on 
the maser morphology and polarization are presented in § 
The results are discussed in § |4] where intrinsic properties of 
the masing regions are derived. This is followed by a summary 
and conclusions in § [5] and § [6] The analysis method and the 
H2O maser models used are presented in AppendixlAl 



2. Observations 

The observations were performed with the NRAO 1 VLBA on 
October 3 2004. The average beam width is * 0.5 x 0.5 mas 
at the frequency of the 616 - 523 rotational transition of H2O, 
22.235080 GHz. We used 4 baseband filters of 1 MHz width, 
which were overlapped to get a total velocity coverage of 
* 44 km s _1 , covering most of the velocity range of the 
H2O masers around the mean velocity of the H2O masers of 
HW2 V y = -11.7 km s " 1 (T9 6). Similar to the observa- 
tions in Vlemmings et aO §002) (hereafter V02) of circum- 
stellar H2O maser polarization, the data were correlated mul- 
tiple times with a correlator averaging time of 8 sec. The ini- 
tial correlation was performed with modest spectral resolu- 
tion (128 channels; 7.8 kHz= 0.1 km s _1 ), which enabled us 
to generate all 4 polarization combinations (RR, LL, RL and 
LR). Two additional correlator runs were performed with high 
spectral resolution (512 channels; 1.95 kHz= 0.027 km s _1 ), 
which therefore only contained the two polarization combina- 
tions RR and LL, to be able to detect the signature of the H2O 
Zeeman splitting across the entire velocity range. The observa- 
tions on Cepheus A HW2 were interspersed with 15 minute ob- 
servations of the polarization calibrator J2202+4216 (BL Lac). 
Including scans on the phase calibrators (3C345 and 3C454.3) 
the total observation time was 8 hours. 



1 The National Radio Astronomy Observatory (NRAO) is a facility 
of the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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2.1. Calibration 

The data analysis pat h is described in detail in V02. It fol- 
lows the method of Kemball et al. ( 1995) and was performed in 
the Astronomical Image Processing Software package (AIPS). 
The calibration steps were performed on the data-set with 
modest spectral resolution. Delay, phase and bandpass cali- 
bration were performed on 3C345, 3C454.3 and J2202+4216. 
Polarization calibration was performed on the polarization cal- 
ibrator J2202 +42 16 (Fig0. Fringe fitting and self-calibration 
were performed on a strong (~80 Jy beam -1 ) maser feature 
(at Vi sr = -15.72 km s _1 ). The calibration solutions were 
then copied and applied to the high spectral resolution data- 
set. Finally, corrections were made for instrumental feed po- 
larization using a range of frequency channels on the maser 
source, in which the expected frequency averaged linear polar- 
ization is close to zero. In order to make a comparison with 
previous results we have used the AIPS task FRMAP in an 
attempt to determine the position of the reference feature be- 
fore any self-calibration or fringe fitting. Though an exact po- 
sition determination was impossible, we found it to be within 
-25 mas of our pointing position (a(72000) = 22' , 56'"17 ! .977 
and 5(72000) = +62°0l'49".419), which was the brightest 
maser feature of the maser region R4 from G03. 

An initial image cube with low resolution (2048 x 2048 pix- 
els of 1 mas) was created from the modest spectral resolution 
data set using the AIPS task IMAGR. In this cube a search was 
performed for maser features and 4 distinct regions with maser 
emission were detected (further labeled I through IV; shown 
in Fig. |2}. For these fields, typically ~100xl00 mas in size, 
IMAGR was used to create high spatial resolution (1024 x 1024 
pixels of 0.09 mas) Stokes I, Q and U image cubes from the 
modest spectral resolution data set. Stokes I and V cubes for 
the same regions were created from the high spectral resolu- 
tion data set. In the high spectral resolution total intensity chan- 
nel maps, the noise ranges from «15 mJy in the channels with 
weak maser features, to «35 mJy when dominated by dynamic 
range effects in the channels with the strongest maser features. 
In the circular polarization polarization maps the rms noise is 
=sl5 mJy. In the lower resolution Stokes Q and U maps the rms 
noise is a; 10 mJy. 

Unfortunately, we found that in a small range of frequency 
channels where a higher frequency band overlaps the neigh- 
boring lower band, cross-talk between the sub-bands resulted 
in unreliable calibration. Although we were able to image 
the masers in those channels (Vi sr between -12.8 km s _1 and 
-13.5 km s _1 as well as between —1.2 km s _1 and -2.5 km s _1 ) 
they were not included in our polarization analysis as the cali- 
bration accuracy was insufficient. 

To calibrate the polarization angle x - 1/2 x atan(f//0 
of the resulting maps, the polarization calibrator J2202+4216 
was mapped using the full 4 MHz bandwidth. The resulting 
map with polarization vectors is shown in Fig. [2 The polar- 
ization vectors were rotated to match the polarization angle 
of J2202+4216 determined in the VLBA polarization calibra- 
tion observations 2 . As our observations were made exactly be- 

2 http://www.aoc. nrao.edu/~ smyers/calibration/ 



tween two of the calibration observations on September 19 and 
October 17 2004 where the polarization angle of J2202+4216 
changed from 41° to 57°, we use the average of 49°. Thus, we 
estimate our polarization angles to contain a possible system- 
atic error of ~8°. 

2.2. Cepheus A HW2 

We detected 4 distinct regions of H2O maser emission between 
Vi sr = -22.5 and 0.5 km s . We did not detect any of the 
maser features with positive velocity from T98 and TOla to a 
limit of «45 mJy. In Fig.|2]we show a 2.3" X 2.3" area around 
HW2 in which the fields where H2O maser emission was de- 
tected are marked. We also indicate the continuum source HW2 
(T96) and the location of previously detected H2O maser not 
visible in our observations. Add itionally, the location of OH 
masers {B arfkiewi cz et all E()05) and 12 GHz methanol masers 
(Min ier et all 1200 ll) are plotted. All offset positions in this pa- 
per are given with respect to reference maser feature position 
at Vis, = -15.72 km s which was earlier found to be within 
25 mas of our pointing center. The accuracy of each individ- 
ual maser feature position can be estimated by Beamsize/SNR, 
which is typically better than ~0.005 mas. In our polariza- 
tion analysis we only considered maser features with intensities 
> 1 Jy. 

3. Results 

3.1. Distribution of the Maser Features 

In Fig.|3]we show the 4 fields in which maser features stronger 
than 1 Jy beam 1 were identified. The hexagonal symbols de- 
noting the maser features are scaled logarithmically by their 
flux density level. We identified 54 maser features although 14 
of those had a velocity located in the ranges that suffered from 
interference as discussed above. The maser features are listed 
in Table^with their positional off-set from the reference maser 
position, peak flux density, radial velocity Vi sr and full width 
half maximum (FWHM) Avi. The positions were determined 
in the frequency channel containing the peak Stokes I emission 
using the AIPS task JMFIT. The masers in Field I, III and IV 
were seen previously (TO lb) while Field II contains a newly 
detected linear maser structure approximate 1" East of HW2 
(assuming a distance of 725 pc the masers are located ~690 AU 
from HW2). The masers in Field I are identified as the masers 
seen in the region labelled R4 of TOlb observed with similar 
lsr velocity. This are also the masers that were hypothesized to 
belong to a rotating disk in G03 and are found over a large ve- 
locity range ( Vi sr between -20 and km s _1 ). As seen in Fig. [2 
the maser structure in Field II is located close to the bright est of 
the 12 GHz methanol maser features detected by M inier etafl 
(2001). However, the methanol masers at Vi sl - = -4.2 km s _1 
are significantly red-shifted with respect to the H2O maser 
structure, which has an average Vi sl - ~ -13.7 km s -1 . The 2 
maser features in Field III at V] sr w -8.5 km s _1 , correspond to 
a small part of the extended maser arc in R5 of TOlb. This arc 
was identified to belong to a spherical shell around an embed- 
ded YSO (Curiel et al., 2002). The fairly weak masers of Field 
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Fig. 2. The Cepheus A HW2 region with related maser features. The crosses indicate the H2O maser features from our obser- 
vations with the boxes labeled I through IV the fields in which the masers where detected. The other H2O maser positions are 
from the V LA and VLBA obs ervations in T96 and T01 respectively. The solid dots are the positions of the 12 GHz me thanol 
maser s from Minier et alJj200ll) (M01) The solid squares are the 1665 and 1667 MHz OH masers observed by Bartki ewicz et alJ 
l l2005l) (B05). The ellipse denotes the position and shape of the HW2 continuum emission at 1.3 cm (T96). 



IV at Vi sr m — 21 km s are located closest to HW2 and likely 
correspond to a few isolated features detected in T96. We did 
not detect any of the masers from the arc-like structures in Rl, 
R2 and R3 of TOlb. The total extent of the region in which we 
detected maser emission is ~ 950 x 790 mas, corresponding to 
690 x 575 AU. 



3.2. Circular Polarization 

Circular polarization between 0.018-2.31% was detected in 14 
of the 40 maser features that did not suffer from the frequency 
band overlap interference. Features that were not analyzed due 
to the interference are marked in Tabled This table also shows 
the circular polarization fraction Py as well as the magnetic 
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Fig. 3. A close-up view of the 4 fields in which we detected H2O maser features. The octagonal symbols are the identified maser 
features scaled logarithmically according to their peak flux density. The maser velocity is indicated by color, note that the color 
scale is different for the 4 fields. A 10 Jy beam 1 symbol is plotted for illustration in the lower left corner of Field IV. The linear 
polarization vectors, scaled logarithmically according to polarization fraction Pi, are over-plotted. For the maser features where 
the Zeeman splitting was detected the magnetic field strength is indicated in mG. 



field strengths along the line of sight with l<x errors or 3<x 
upper limits determined by comparing the line width and cir- 
cular polarization with models of non-LTE radiative transfer 
in the magnetized H2O molecules (Appendix IaV As the lcr 
errors include both the formal fitting uncertainties as well as 
the contribution of the error in the model Av t h (thermal line 
width) and T^AQ. (emerging maser brightness temperature in 
K sr), the magnetic field strength can occasionally be < 3<x, 
even though the circular polarization signal has a SNR higher 
than 3. The table also includes the best fit model values for 
Av t h and 7t,A£2, where the emerging brightness temperature has 
been scaled with maser decay and cross-relaxation rate as de- 
scribed in Appendix|A| The errors on these are estimated there 
to be 0.3 km s _1 in Av t h and 0.4 on log(7bAf2). As the lack of 
circular polarization introduces an additional free parameter in 



the model fitting, significantly increasing the Av t h and T^AO, er- 
rors, we do not fit for maser features that do not show circular 
polarization. The magnetic field strength ranges from several 
tens of mG in Fields II and III to several hundred mG in Field 
I and is seen to switch direction on small scales in both Field 
I and II. Note that a positive magnetic field values indicates a 
field pointing away from the the observer. Total intensity (I) 
and circular polarization (V) spectra of several of the maser 
features are shown in Figs. |4] and |5] The spectra include the 
best fit model for the circular polarization. 

3.3. Linear Polarization 

In addition to the circular polarization, we detected linear po- 
larization in approximately 50% of our maser features. The 
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Fig. 4. Total power (I) and V-spectra for selected maser features of Field I. Additionally, the linear polarized flux density, V(6 2 + 
U 2 ), is shown when detected. The flux densities are given in Jy beam 1 . The thick solid line in the bottom panel shows the best 
non-LTE model fit to the circular polarization V. The V-spectrum is adjusted by removing a scaled down version of the total 
power spectrum as indicated in AppendixlAl 
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Fig. 5. Similar to Fig.|4]for selected H2O masers in Field II and both features in Field III. 



fractional linear polarization P/ is given in Table [2 Figs. 0] 
and[5]also show several linear polarization spectra. Table[^lists 
the weighted mean polarization vector position angle (x) deter- 
mined over the maser FWHM for the linearly polarized maser 
features with corresponding rms error. The weights are deter- 
mined using the formal errors on^ due to thermal noise, which 



are g iven by <r x - 0.5 <t p /P x 180° /n (Ward le & Kronbergl 
119741) . Here P and cr P are the polarization intensity and cor- 
responding rms error respectively. Fig. [3] shows the linear po- 
larization vectors scaled logarithmically according to fractional 
polarization. 



Table 1. Results 
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Feature 


RA 

offset 


Dec 
offset 


Peak Flux 
Density (I) 




Av L 


P 


(x) 


Pv 


By" 


Av th " 


io g (r b An) a 




(mas) 


(mas) 


(Jy beam ) 


km s -1 


km s 


(%) 


o 


(xl0~ 3 ) 


(mG) 


km s 




La* 


44.213 


-57.325 


10.32 


-1.86 


0.53 


- 


- 


_ 


_ 


- 


- 


Lb 


43.289 


-56.585 


16.79 


-3.83 


0.59 


0.28 ±0.12 


13 + 8 


2.6 


62+ 12 


1.8 


10.5 


I.c 


42.750 


-56.510 


9.83 


-4.12 


0.70 


0.44 ± 0.06 


35 + 8 


9.5 


-290 ± 47 


2.0 


10.7 


I.d" 


42.290 


-47.644 


19.09 


-4.28 


0.58 


0.78 ± 0.06 


-46 + 3 


12.7 


-279 ± 69 






Le 


42.229 


-48.655 


75.18 


-3.96 


0.78 


0.64 ± 0.04 


-47 + 8 


10.3 


205 ± 40 


1.0 


10.0 


I.P 


34.428 


1.488 


33.35 


-13.04 


0.48 


- 


- 


_ 


_ 


- 




!g* 


33.811 


1.977 


27.49 


-12.97 


0.53 


- 


- 






- 




I.h* 


32.322 


-46.771 


87.44 


-1.94 


0.62 


_ 


_ 


_ 


_ 


_ 


_ 


Li* 
I-j* 


31.031 
29.643 


-47.121 
-47.557 


12.81 
8.49 


-1.68 
-1.20 


1.01 
0.94 


- 
_ 


- 
_ 


_ 


_ 


- 
_ 


_ 


Lk 


1.963 


-1.027 


3.48 


-18.83 


0.57 


< 0.86 


- 


- 


< 206 


- 


- 


LI 


1.040 


-2.099 


3.04 


-19.04 


0.51 


< 0.99 


- 


23.1 


527 ± 109 


2.0 


9.9 


Lm 


0.082 


-3.303 


11.47 


-19.44 


0.50 


< 0.26 


- 


6.8 


135 ± 26 


1.7 


10.3 


Ln 


0.000 


0.000 


78.94 


-15.72 


0.52 


< 0.04 


- 


7.2 


148 ± 34 


1.7 


10.4 


1.0 


-0.820 


-4.152 


4.02 


-19.28 


0.59 


1.17 + 0.19 


8 + 9 


- 


< 256 


- 


- 


Lp 


-0.871 


-0.183 


61.63 


-15.89 


0.51 


0.42 ± 0.02 


58 + 2 


6.9 


-150 + 42 


2.0 


9.8 


I.q** 


-1.486 


-1.359 


46.82 


-16.78 


0.68 


0.59 ± 0.26 


39 + 2 


6.8 


-203 ± 71 


- 


- 


E.a 


925.674 


13.463 


18.07 


-14.91 


0.51 


0.33 ± 0.09 


-41+5 


_ 


< 22 


- 




n.b 


924.445 


13.075 


44.81 


-14.93 


0.49 


< 0.07 


- 


3.5 


-54 + 9 


1.5 


9.8 


n.c 


906.332 


6.070 


2.76 


-14.25 


0.48 


< 1.09 


- 




< 140 


- 




n.d 


896.763 


3.990 


15.24 


-14.06 


0.51 


0.52 + 0.01 


9 + 7 




< 34 


- 




n.e 


894.765 


3.352 


21.08 


-13.99 


0.59 


0.37 ± 0.05 


-12+ 10 


_ 


< 28 


_ 


_ 


n.f 


892.838 


2.587 


24.33 


-14.09 


0.46 


0.21 ± 0.05 


38+ 12 




< 30 


- 




n.g 


890.204 


1.947 


3.79 


-13.99 


0.33 


< 0.79 


- 




< 88 


- 




n.h 


888.068 


3.166 


3.20 


-13.55 


0.40 


3.44 ± 1.13 


-78 + 4 


_ 


< 125 


_ 


_ 


Hi 


887.387 


1.437 


33.64 


-13.70 


0.50 


0.15 + 0.02 


-50+ 16 




< 15 


- 




n.j 


885.188 


0.222 


8.03 


-13.60 


0.47 


1.05 + 0.10 


6 + 4 




< 59 


- 




n.k* 


882.617 


-0.579 


11.25 


-13.47 


0.53 


_ 


_ 


_ 


_ 


_ 


_ 


n.i* 


881.495 


-0.880 


6.95 


-13.31 


0.47 


- 


- 






- 




Il.m* 


880.610 


0.001 


5.86 


-13.26 


0.46 


- 


- 






- 




E.n* 


876.069 


-4.921 


5.25 


-12.81 


0.40 


_ 


_ 


_ 


_ 


_ 


_ 


n.o* 


875.052 


-5.871 


5.03 


-12.81 


0.44 


- 


- 




_ 


- 




n. P 


870.917 


-11.180 


2.48 


-12.57 


0.47 


< 1.21 


- 


- 


< 190 


- 


- 


n. q 


869.659 


-11.433 


5.84 


-12.65 


0.50 


< 0.51 


- 


3.5 


66 + 33 


1.8 


10.2 


Il.r* 


858.940 


-15.758 


7.85 


-12.91 


0.41 


- 


- 


- 


- 


- 


- 


II.s* 


857.685 


-16.246 


8.92 


-12.91 


0.40 


- 


- 


- 


- 


- 


- 


n.t 


855.098 


-11.502 


5.21 


-14.51 


0.55 


0.41 ± 0.07 


86+ 19 


- 


< 127 


- 


- 


n.u* 


854.380 


-17.668 


56.09 


-12.83 


0.42 


- 


- 


- 


- 


- 


- 


n.v 


853.907 


-18.350 


23.43 


-12.78 


0.41 


0.34 ± 0.09 


84 + 7 


5.7 


69+ 11 


1.1 


9.7 


III. a 


143.668 


-530.929 


21.04 


-8.25 


0.58 


10.8 + 0.9 


66+1 


1.8 


33+10 


1.2 


10.7 


m.b 


142.096 


-532.309 


6.02 


-8.52 


0.57 


5.0 + 0.8 


62 + 2 


6.9 


128 ± 36 


1.3 


10.5 


IV. a 


60.131 


229.423 


1.49 


-21.28 


0.62 


< 2.01 






< 520 






IV.b 


60.025 


226.307 


2.28 


-21.07 


0.65 


1.35 ± 0.24 


64 + 2 




< 356 






IV.c 


58.833 


240.075 


1.38 


-20.75 


0.86 


< 2.17 






< 779 







" Best fit results for the magnetic field strength Etjj along the line of sight (mG), intrinsic maser thermal width Av lh (km s~') 
and emerging brightness temperature Tt,AQ (K sr) derived as described in AppendixlAl 
* suffer from interference (see text) 
** no direct fit possible (see AppendixlAl 



The strongest linear polarization (~11%) was detected 
on the brightest maser feature in Field III, but on average 
Pi ~0.5%. In Fig. [6] we show a channel map of the 2 maser 
features detected in Field III including their polarization vec- 
tors. We do not find any relation between maser brightness and 
fractional linear polarization. 



4. Discussion 



Before discussing the polarization results we first determine 
several intrinsic properties of the masers that are needed for 
the further analysis of the linear and circular polarization. We 
also discuss the maser morphology in the 4 fields. 
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Fig. 6. Channel maps of linear polarization of the elongated H2O maser feature of Field III which has the highest linear polariza- 
tion fraction. The bars show the strength and orientation of the polarization vectors. 



4.1. Intrinsic Thermal Width, Brightness Temperatures 
and Maser Beaming 

As the model results give the intrinsic thermal width Av t h in 
the maser region, we can use it to estimate the temperature. 
Although the error on Av t h i s relatively large due to veloc- 
ity gr adients along the maser (Vlemmings & van La ngevelde. 
2005), we find that on average, Avth, and correspondingly the 
temperature, is greater in Field I than in the outlying fields II 
and III. While in Field I T ~ 1 150 K, in Field II and III the cor- 
responding temperature is closer to 750 K. These temperatures 
are an indication that the masers originate in a C-type (non- 
dissociative) shock instead of a J-type (dissociative) shock. In 
the latter, the H2O masers have been found to originate in a rel- 
atively narrow range of temperatures near 400 K (with 500 K 
as a conservative upper bound) at which hydrogen molecules 
recombine. In contrast, in C-type shocks, the H2O masers can 
occur in gas with temperatures up to ~3000 K provided th e 
shock velocity v., > 10 km s~' ( Kaufman & Neufeldt 1 1996). 
The C-shock origin of the masers in Field I is in agreement 
with the model in G03 where the masers originate in a C-shock 
expanding though a circumstellar disk. 

In addition to Av t h the models also provide an estimate of 
the emerging brightness temperature 7bAQ. This can be com- 
pared with the values determined from the measurements of 



the maser flux density and feature sizes. We find that in Field 
I the majority of the maser features are unresolved. Taking 
0.4 mas as the typical size of a H2O maser feature, we de- 
rive a the brightness temperature of 7b « 1.4x10" K for a 
feature of 10 Jy beam -1 . Thus, our strongest maser feature 
in Field I has 7b ^l.lxlO 12 K. In the other fields, several of 
the masers are marginally resolved, with typical feature sizes 
of ~0.6 mas, corresponding to ~7.5xl0 12 cm. This implies, 
for the strongest 54 Jy beam -1 maser feature in those regions, 
7b ~3.4xlO n K. Comparing these values with the emerging 
brightness temperatures T^AD. from our models yields an esti- 
mate for the beaming solid angle AO. In Field I, with an av- 
erage (TbAQ) »1.8xl0 1() we find, for the maser features with 
circular polarization, AQ »7xl0 -3 -3xl0 _1 sr. However, as the 
features are unresolved the beaming angle may be overesti- 
mated. The masers in Field II show a similar range of beam- 
ing angle, with AO a;2xl0 -2 -4xl0 -1 sr while the beaming of 
the maser in Field III is much less pronounced, as AQ »0.5. 
In a tubular geometry the maser beaming AQ. m (d/l) 2 , where 
d and / are the transverse size and length of the tube respec- 
tively, this implies that, assuming d is approximately the size 
of the maser features, the maser lengths are ~l-6xl0 13 cm. In 
Field III, the beaming angle is similar to what is expec ted for a 
spherical maser that is approaching saturation ( Eli tzuii Tl 994l) . 



W.H.T. Vlemmings et al.: Magnetic Fields in Cepheus A 



9 




Fig. 7. a) The masers of Field I with the disk-model of G03. 
We fitted a Right Ascension and Declination offset of the disk 
center (denoted by the plus sign) as well as a Radius {R4 = 
34 ± 4 mas). The inclination angle of 50° and P.A. of 142° 
were taken from G03. The diagonal cross is the disk center 
position determined by G03 which has an error of ^27 mas 
in each coordinate. The vectors on the maser features are the 
polarization vectors, which for most of the features is expected 
to be parallel to the magnetic field direction (see § !4.3t . b) The 
masers of Field III with the expanding shell model of TO la. 
The solid line is the shell using the updated proper motion and 
expansion velocity parameters of G03. The dashed lines are 
the 3cr confidence interval. The vector on the maser features 
indicate the magnetic field direction (see § I4.3> . 

We now compare our measured and derived maser bright- 
ness temperatures with the maser brightness temperature 7^ 
at the onset of saturation when the ratio between maser rate 
of stimulated emission (R m ) and t he maser decay rate (F), 
R m /T ~1. Using the expression from lReid & Moranl 1^988): 

T S AQ = hvT4nl2k B A, (1) 

where h is the Planck constant and kg the Boltzmann constant, 
v is the maser frequency and A — 2 x 1 0~ 9 s' 1 is the 22 G Hz 
H?Q maser spontaneous emission rate JGoldreich & Keelevl 
I1972I) . For F = 1 s 1 we thus find T s ACl = 3.4 X 10 9 K sr. 
Nearly full saturation is reached when R m /T ~100, for T\,AQ. « 
3.4 X 10 11 K sr. This indicates that in Field II the masers are 
likely mostly unsaturated, while those in Field I are in the onset 
of the saturation regime. In Field III the masers are almost fully 
saturated, which is consistent with their strong linear polariza- 
tion (see Appendix lAl below). When saturated, the maser ra- 
diative transfer equation can be approximated by 7b /7s ~ gol 
where go is the maser gain at line center for the unsaturated 
regime. For the masers in Field III we then find gol ^8, indicat- 
ing that for / «1.5xl0 13 cm estimated from the beaming angle, 
go *5xl0 -13 cm" 1 . 

4.2. H 2 Maser Morphology 

As can be seen in Figs. |2 and [3] the H2O masers around 
Cepheus A HW2 show a large variety of structures. In our ob- 
servations several of the maser structures found in TOlb and 
G03 were not detected, even though our sensitivity is within 
a factor of 2 of those of TOlb (~6 mJy) and G03 (~25 mJy). 



Fig. 8. The velocity of the H2O maser features making up the 
filament detected in Field II vs. angular off-set from the cen- 
ter of the filament. The solid and open circles are the maser 
features that have detected linear polarization while the open 
squares are the features for which we determined upper limits 
to Pi. The crosses are the features that are affected by the in- 
terference described in §|2]and which were excluded from our 
polarization analysis. The open circle is a feature that likely 
does not belong to the filament and which has been excluded 
in the subsequent analysis. The solid line is a best fit relation 
between the maser velocity and position off-set. 



We did detected a strong linear maser structure in Field II that 
was not observed in the previous observations. The changes 
in observed morphology are likely due to the rapid variability 
of the H2O masers of Cepheus A, which show variat i ons o n 
timescales as short as a few days (Rowland & Cohen, 1986). 
Here we discuss the masers of the 4 distinct regions detected in 
our observations. 

Field I: The H2O masers in this field show the most com- 
plex structure. The masers are located ^150 AU on the sky 
south of the continuum source HW2 and have been previ- 
ously detected in T96, TOlb and G03. In TOlb this region was 
named R4 and it was proposed that the masers of the NW 
corner (named R4-A) originate in a bow-shock structure pro- 
duced by the wind of an undetected protostar near R4-A. The 
features in the SE could be connected to R4-A and produced 
by a shock moving ~4-7 km s _1 to the NE. In G03, where 
these masers were observed with MERLIN, the masers are, in- 
stead of in a bowshock, hypothesized to occur in an expand- 
ing Shockwave in a rotating proto-stellar disk enclosing a cen- 
tral mass of ~3 M©. Here we only detect several bright fea- 
tures making up an incomplete disk with Vi sr w -4.0, - 1 .5 and 
-19.0 km s . We fitted our maser feature positions to the disk 
proposed by G03 using a flux density weighted least square 
method. Keeping the inclination and position angle fixed with 
the G03 model values (50° and 142° respectively) a fit was 
made for the Right Ascension and Declination offset of the 
disk center and for the disk Radius (Rd)- The result is shown in 
Fig- Et- Our disk center offset position is only ~1 mas SW of 
the position determined by G03 while the error on the reference 
position determination in this paper combined with that of G03 
is estimated to be ~27 mas. Our fitted disk radius Rd = 34 mas. 
Considering we only detect a small part of the disk and since 
the masers in the SE corner are spread over a large area we es- 
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timate the systematic error in our radius determination to be 
4 mas, larger than the formal fitting uncertainty of ~1 mas. 
Comparing with the radius determined at Epoch 2000.27 
by G03 (R d = 38.1 ± 0.1 mas) we find that the disk has not ex- 
panded in the 4 years after their observations. It possibly even 
decreased in radius. In G03 it was found that the expansion 
velocity decreased from 30-40 km s" 1 in 1996 to ~13 km s" 1 
in 2000. This strong deceleration apparently has continued and 
may be due to mass loading of the disk as matter is swept up 
during the expanding Shockwave. As a result a stationary shock 
may have formed where the circumstellar outflow collides with 
the much denser surrounding medium. This could also explain 
the disappearance of the brightest disk masers observed in G03, 
since for higher shock number densities (> 10 12 cm" 3 ) the 
masers will be quenched. 

Field II: The H2O masers in Field II make up a newly dis- 
covered filamentary structure ~690 AU East of HW2 at a posi- 
tion angle (PA) of 66.0° ± 0.2° and with a length of -60 AU. 
This structure also nearly coincides with 12 GHz methanol 
masers (at different Vi sl ) located ~40+10 mas to the NW, which 
show a linear structure with similar PA. As seen in Fig.[8]there 
is a velocity gradient along the filament from — 15 km s" 1 in 
the NE to — 12.5 km s" 1 in the SW. The maser structure bears 
resemblance to the masers in Rl, R2 and R3 of T01 found to- 
wards the West of HW2 with a similar PA, although the masers 
in Field II are all blue-shifted with respect to the systematic ve- 
locity of HW2 while those in Rl, R2 and R3 were red-shifted. 
The masers are too far East to be considered part of the rotating 
maser disk around HW2 which is thought to have a radius of 
300 AU (T96). The elongated appearance of the Field II maser 
structure sugg ests a shock ed nature as expected from maser 
theory felitzur et allll989h . Although it is located at a signifi- 
cant distance from HW2 we suggest that the maser structure is 
due to the interaction of the HW2 outflow with the circumstel- 
lar molecular cloud medium. Then, similarly as for the Western 
Rl features in TOlb, the velocity shift of ~2.5 km s" 1 along the 
maser filament could be due to acceleration of maser gas by 
the YSO outflow. If the masers are indeed created by shocks in- 
duced by the HW2 outflow this would indicate that at ~690 AU 
the outflow has an opening angle of ~1 15°, similar to the open- 
ing angle of ~110° estimated for the Rl masers at 150 AU in 
T01. 

Field III: The 2 maser features detected in Field III, approx- 
imately 550 AU South of HW2 at V lsr *-8.5km s"\ likely be- 
long to the arc structure R5 described in TOla and TOlb. These 
masers are thought to be part of a spherical she ll surround- 
ing a p rotostar that has possibly been identified in lCuriel et alJ 
(2002). We do not detect the long maser arc seen in TOla and 
TOlb. While the brightest maser feature in our observations 
of Field III is ^20 Jy beam" 1 , the brightest maser features of 
R5 in TOla and TOlb (separated by 8.5 yr with respect to our 
observations) had flux densities of ^200 Jy beam" 1 . However, 
the PA (~41°) of the extended maser structure of our obser- 
vations (seen in Fig.[6ji agrees with the direction of the maser 
curve. Fig. 05 shows the maser shell from TOla extrapolated 
in time using the updated shell parameters determined in G03. 
The curves indicate a shell expansion of 2.5 ±0.1 mas yr" 1 
and a motion of the expansion center of 1.4 ± 0.1 mas yr" 1 



toward PA 126°. The near-perfect alignment of the expanding 
shell with the maser features is remarkable, as we earlier esti- 
mated our reference position to be accurate to ~25 mas. This 
likely indicates that we underestimated our positional accuracy 
and that the actual accuracy is better than ~10 mas. In addi- 
tion, our results indicate that the maser shell has been freely 
expanding during the past 8.5 years without any indication of 
deceleration. 

Field IV: The masers in Field IV are weak and are aligned 
at a PA ~ -6°. They are located within -75 AU of HW2 at 
Vlsr »-2 1.0km s" 1 and are probably part of the rotating maser 
disk around HW2 proposed in T96. 

4.3. Linear Polarization 

Linear polarization is often affected by Faraday Rotation due 
to free electrons along the line of sight through the interstel- 
lar medium. However, the Faraday rotation induced in a typ- 
ical molecular cloud with fairly strong magnetic field (size 
D ~Q.l pc, electron density n e ~1 cm" 3 and B\\ ~1 mG) is only 
~ 0.9° at 22.235 GHz. The rotation induced in the extreme con- 
dition of a highly magnetized maser cloud (up to 1 G) is similar 
or less. As no compact HII regions, in which Faraday rotation 
could be significant, are located in front of the maser features, 
we can safely assume the measured ;f is not affected by Faraday 
rotation. 

As discussed in Appendix|X| the polarization vectors deter- 
mined from polarization observations of masers in a magnetic 
field are either parallel or perpendicular to the magnetic field 
lines. Thus, the polarization vectors contain information on the 
morphology of the magnetic field but suffer from a 90° degen- 
eracy. The fractional linear polarization depends on the maser 
saturation level as well as the magnetic field angle 6. Thus, we 
can use the measurements of Pi together with our model re- 
sults for the saturation level (through the emerging brightness 
temperatures) to lift the degeneracy between the polarization 
vectors and the direction of the magnetic field for several of 
our maser features. 

The polarization vectors observed in the maser fields 
around Cepheus A HW2 are shown in Fig.[3]while p t is listed in 
Tabled The strongest linear polarization Pi » 1 1 % was found 
in Field III. This is consistent with the fact that the brightness 
temperature analysis concluded that the masers in this field are 
saturated. Using the brightness temperature determined from 
the models, adjusted for the difference in maser decay and 
cross-relaxation rate as described in Appendix|A] we find using 
Fig. IO that for the masers i n Field III, 65° < 6 < 70°. As was 
shown in IVlemmin gs (2005, hereafter V05) this is the mag- 
netic field angle in the unsaturated (or least saturated) maser 
core. Thus, as 6 > 6 cl -i t , the magnetic field direction is perpen- 
dicular to the polarization vectors. As can be seen in Fig. 0> this 
means the magnetic field, at PA~155°, is perpendicular to the 
expanding shell found in TOla and thus radial from the central 
embedded proto-star. 

In Field IV fairly strong linear polarization was detected in 
one of the weak maser features. As the brightness temperature 
of these masers is relatively low and they are unlikely to be 



W.H.T. Vlemmings et al.: Magnetic Fields in Cepheus A 



11 




Fig. 9. (Left) x f° r the H2O maser filament features in Field II that have measured linear polarization vs. their velocity. The 
connected dots belong to the individual features which are labeled with the feature identifier. We excluded the seemingly unrelated 
feature denoted by the open circle in Fig|8] The thick solid lines (with a 180° ambiguity) is a linear fit to the change in^- along 
the feature. The thick dashed lines indicates the linear fit when including the 90° flip in x with respect to the magnetic field 
direction when the magnetic field angle to the line-of-sight 9 becomes larger than 9 CT n ~55°. For the features connected with the 
solid lines we expect the x to be parallel to the direction of the magnetic field while for those connected with dashed lines, x is 
perpendicular to the magnetic field. (Right) x for the H2O maser in Field I, III and IV that have measured linear polarization vs. 
their velocity. The connected dots belong to the individual features that are labeled with their corresponding identifier. The boxes 
are labeled with the field number. 



saturated, Pi = 1.35% indicates that the magnetic field angle 
9 > 70°. Thus also in Field IV the magnetic field direction is 
perpendicular to the polarization vector with a PA~154°, more 
or less along the large scale maser disk proposed in T96 with a 
PA~135° and radial toward HW2. 

The H2O masers in the circumstellar disk of Field I are 
found to have Pi < 1 %, consistent with their being only slightly 
saturated. With 7b determined earlier, we find, again using 
Fig. IA.ll that 9 is either close to 6* CI i t or 9 < 25°. As seen in 
Fig-Et, the polarization vectors mostly lie along the disk cur- 
vature for most features except I.c and I.d, where we could be 
seeing a 90° flip. If we assume that for all features, except I.c 
and I.d, 9 < 6* cl ; t , the magnetic field in Field I lies along the 
disk. However, if most features have 9 > # cr jt except for I.c 
and I.d, the magnetic field is radial in the H2O maser region of 
the rotating disk. Additionally, as in either case, 9 is close to 
9 a h = 55° the magnetic field angle with respect to the line-of- 
sight is similar to the disk inclination axis, which would imply 
the magnetic field lies in the plane of the disk. However, as 
seen in Table [I] the magnetic field direction changes between 
the neighboring maser features making a large scale alignment 
unlikely. 

Now we show that the polarization characteristics of the 
masers in Field II are consistent with the interaction between a 



radial magnetic field in the outflow of HW2 and a magnetic 
field perpendicular to the Galactic plane in the surrounding 
molecular cloud. The fractional polarization of the maser in 
the filamentary structure of Field II is on average slightly less 
than that in Field I. This is expected since the masers in Field II 
were found to be unsaturated. The high polarization of feature 
Il.h likely indicates that there 9 is close to 90°. As seen in Fig.[3J 
there is evidence of a gradient in polarization angles along the 
maser filament. In addition to the gradient along the maser fil- 
ament, we see in the left panel of Fig. [9] that the polarization 
angle x rotates across individ ual maser features, sim ilar to that 
seen in the cocoon masers of iLeppanen et alJ lll998l) . Such ro- 
tation of x is not observed for any of the maser features of the 
other fields shown in the right panels of Fig. [5] The variation of 
X with velocity can be described with a linear gradient, using 
an flux density weighted least square method allowing for the 
90° flip inx that occurs when 9 > cl ; t . We find that;^ increases 

linearly from 50° on the maser feature in the NE to ~90° 

on the feature in the SW. This implies that Il.h and II. i undergo 
the 90° flip which was already expected for Il.h due to its high 
polarization. 

Similar to the model for the variation in polarization an- 
gle x, we have constructed a model for the variation of the an- 
gle between the magnetic field direction and maser propagation 
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axis 8 along the maser filament. The model, shown in Fig. 1101 
is fully consistent with the maser brightness temperatures and 
the fractional linear polarizations as well as the inferred 90° 
flip of polarization angle. We have determined 8 and its error 
bars from the relation between P/ and 8 for unsaturated masers 
shown in Fig. IA.ll (Note that to the accomodate the direction 
change of the magnetic field, the model 8 ranges from 0° to 
180°, with the direction change occurring at 8 = 90°). We find 
that the polarization measurements are consistent with an initial 
slow change in 8 until halfway along the filament the magnetic 
field changes sign over ^10 mas. 

Combining the x an d 8 variation models we thus find that 
at the NE corner of the filament the magnetic field is pointing 
toward us with 8 -10° and PA 50°, while in the SW cor- 
ner it is pointing away from us with 8 ~2Q° and PA~90°. We 
interpret this change of the magnetic field as being due to the 
interaction of the magnetic field related to the HW2 YSO out- 
flow and that related to the surrounding medium. Although the 
magnetic field str ucture in the Cepheus A region is complex, 

1 Jones et all (|2004) find using Near Infrared imaging polarime- 
try, that the large scale field threading Cepheus A is almost 
perpendicular to the galactic plane with PA= 125°. This cor- 
responds to PA=-55°, consistent wit h y - —50° fo und in the 
NE corner of the maser filament. |jones et al] J2004h also argue 
that the magnetic field in the HW2 outflow is radial with re- 
spect to HW2. This implies an angle of 1 15° at the location of 
Field II, very consistent with the PA in the SW corner of the 
H2O maser filament, especially as we are probably not probing 
the full polarization vector rotation along the filament. 

4.4. The Magnetic Field 

AAA. The Magnetic Field Strength 

The magnetic field strength was determined from circular po- 
larization measurements for 14 maser features in Field I, II and 
III. In Field I the magnetic field strength varies from B = -290 
to 527 mG, while in the other fields we find B between -54 
and 128 mG. While these latter magnetic field strengths are 
comparable to the typical field strength determined from other 
H2O maser observations (10-100 mG), the field strengths de- 
termined in Field I are several times higher. However, previous 
observations wer e typic ally performed using single-dish (e.g. 
iFiebig & Giisterl Pi 989 i or lower resolution interferometers 
(e.g lSarma*etaiIl2002h and due to blending of a large num- 
ber of maser features, the magnetic field strength determined 
with single-dish and VLA observations could be more than a 
factor of 2 smaller than the actual field strengths JSarma et all 
l200lb . The only other H2O mas er magnetic field st rength for 
Cepheus A was determined by Sar ma et al] J2002I) with the 
VLA, who found B =-3.2 mG for a feature located more than 

2 arcsec East of our observed maser region. 

An additional complication to the accurate determination 
of magnetic field strengths is the occurrence of velocity gra- 
dients along the maser amplification path. This was investi- 
gated in V05, where it was found that for velocity gradients 
of ~1.5 km s _1 along the maser, the magnetic field could be 
underestimated by ~40%. From a total intensity line profile 



1 1 ' ' l_ 






" f i 


; — 

VL 


■ 1 1 1 1 1 1 






1 1 1 1 1 




40 20 -20 -40 
Off-set (mas) 


1 , 1 1 1 , , 1 1 1 1 , 1 




Fig. 10. (bottom) The magnetic field angle 8 estimated from 
the fractional polarization P/ measurement. The vertical dashed 
line indicates 6> cl i t . The solid dots are the measured Pi of the 
Field II maser features and the open squares are upper limits. 
The error bars in 8 are determined from the allowed range of 9 
in the Pi vs. 8 models of V05 shown for masers with emerging 
brightness temperatures T\,AQ. = IQ 9 and 10 10 K sr. (lower and 
upper solid line respectively), (top) A model for the change 
of magnetic field angle 8 along the H2O maser filament. The 
thick solid line is the proposed model. The solid horizontal line 
indicates where the magnetic field direction changes. Between 
the dashed horizontal lines, which denote 8 cl -i t , the polarization 
vectors are perpendicular to the direction of the magnetic field. 
Above and below 8 CI i t the polarization vectors are parallel to 
the magnetic field. 



analysis as described in Vlemmings & van Langeveldel J2005) 
we estimate the typical velocity gradient for our masers to be 
~1 km s . For partly or fully saturated masers with Av t h > 
1.5 km s _1 , V05 finds that the magnetic field is overestimated 
by not more than a few percent. However, for the unsaturated 
masers in Field II the field strengths have most likely been un- 
derestimated by ~30%. 

The magnetic field dependence on 8 introduces further un- 
certainties. While for low brightness temperature masers B is 
straightforwardly dependent on cos 8, this relation breaks down 
fo r higher brightness temperat ures. This was first investigated 
in iNedoluh a & Wats onl d 19921) and later shown in more de- 
tail in Watson & Wvld (2001) for masing involving low angu- 
lar momentum transitions. The specific case for the 22 GHz 
H2O masers was again shown in V02 and their figure 7 is re- 
produced here as Fig. IA.2l In the figure we see that for increas- 
ing saturation there is a large range of 8 where the magnetic 
field is actually overestimated. As we have been able to esti- 
mate 8 for several of the observed masers we can also estimate 
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the influence on the magnetic field strength. For the masers in 
Field I 8 is thought to be close to fl cr i t , while the masers are sat- 
urating. This means the actual magnetic field is approximately 
10% higher than the measured field strength. Thus we estimate 
the average field in Field I to be ~250 mG with a maximum 
of ~650 mG. The masers in Field II are however mostly unsat- 
urated and as a result \B\ = By/ cosfl. As we only detected a 
magnetic field strength at the edges of the filament, where we 
estimated 6 to be between 10° and 20°, the field in the NE of the 
filament is ~55 mG pointing toward us while it is ~70 mG and 
pointing away from us in the SW. Finally, the masers in Field 
III were found to have 65° < 6 < 70°. As they are saturated 
the magnetic field strength is likely ~20% less than determined 
from our fits, indicating that |B| = 30-100 mG. 

Aside from the large magnetic field strength, the maser 
structure in Field I is also characterized by field reversals on 
small scales. The magnetic field is found to reverse over less 
than 0.1 mas, which corresponds to ~ 10 12 cm. This argues 
against a large scale alignment of the magnetic field with the 
maser disk. The magnetic field is likely frozen into high den- 
sity maser clumps in a turbulent medium. If the masers exist 
in a shocked region where the magnetic pressure supports the 
cloud and dominates the gas pressure higher magnetic fields 
can b e obtained. Using formula 4.5 from lKaufman & Neu feld 
ill996l) . 

B ~ 80(-/^) 1/2 (— ^— p) mG, (2) 
10 8 cm- 3 10 km s" 1 

where «o is the pre-shock H2 density and v s is the shock ve- 
locity, we find that for a shock velocity v s = 10 km s" 1 as 
estimated for Field I, a magnetic field B - 600 mG can be 
reached if the pre-shock number density «o = 5.6 x 10 9 cm" 1 . 
Estimating the pre-shock m agnetic field using the empirical 
relation of ICrutcheJ {T999) B oc « 47 from the density and 
magnetic field found at the edge of NH^ molecular clouds 
(B = 0.3 mG, n = 2 x 10 4 cm" 3 : iGarav et all 1 19961) yields 
Bq xT00 mG, which is almost 2 orders of magnitude larger 
than the typical pre-shock magnetic field strength (~1 mG). 
Also, when determining the number density of hydro gen in the 
shocke d H2O maser region using the relation from ICrutcheJ 
Jl999l) . the magnetic fields imply densities = 5 X 10 9 - 
2 x 10 11 cm -3 . While the low-end values for n are reason- 
able for H2O masers, the high end (> 10 10 cm" 3 ) is unlikely, 
as such high densities quench the maser population inversion. 
Thus, the magnetic field strength in the pre-shock medium of 
the proto-stellar disk is likely enhanced by the pressence of a 
nearby magnetic dynamo. 

Using Eq.|2]to estimate the pre-shock number density near 
the maser filament in Field II yields, assuming Vs — 13 km s" 1 
similar to the shock velocity in Rl to R3 of TOlb, «o = 3x 10 7 - 
1 x 10 8 cm" 3 . Scaling with B oc n this implies, for the pre- 
shock magnetic field Bo ~10— 15 mG, similar to the magnetic 
field determin ed for co mparable densities in the OH masers 
of Cepheus A (Bartkie wicz et alj.l2005l) . For the number den- 
sity in the shocked region this implies = 3.5 x 10 8 - 
4.7 x 10 9 cm" 3 , typical for H2O masers. 



4.4.2. The influence of the Magnetic Field 

We now examine the influence of the magnetic field on the 
molecular outflow around HW2. When the magnetic field pres- 
sure becomes equal to the dynamic pressure in the outflow 
the magnetic field will be able to influence or even control 
the molecular outflow. Defining B cr \ t the critical magnetic field 
where the dynamic and magnetic pressure are equal, we find 

Bcri, = (87rpv 2 ) 1/2 , (3) 

where p and v are the density and velocity of the maser medium 
respectively. Assuming an outflow velocity of ~13 km s" 1 
we find B cr ; t ^30, 100 and 350 mG for number densities of 
riH 2 = 10 8 , 10 9 and 10 10 cm" 3 respectively. This means that in 
all the H2O maser regions where we measured the magnetic 
field strength the magnetic pressure is approxima tely equal to 
the dy namic pressure, as was previously found in lSarma et alJ 
(2002). As OH maser polarization observations indicate that 
this also holds in the lower density pre-shock regions, we con- 
clude that the magnetic field strength is capable of controlling 
the outflow dynamics. 

5. Summary 

Using polarimetric VLBA observations of the H2O masers 
around Cepheus A HW2 we have been able to measure the 
magnetic field strength and direction in great detail at sub-AU 
scales. We detected H2O masers over an area of ~1 " in 4 dis- 
tinct fields. For each of the fields we derived physical properties 
and several intrinsic properties of the masers. 

Field I: The H2O masers in this field occur in what was pro- 
posed in G03 to be a spherical Shockwave expanding through 
a circumstellar disk. We find that between the G03 MERLIN 
observations in 2000 and our observations the maser ring has 
not expanded and conclude that the expanding Shockwave has 
been severely decelerated, possibly due to mass-loading. From 
our maser models and the measured brightness temperatures 
we find that the typical maser beaming angle in this field is 
~5xl0~ 2 sr implying maser amplification lengths of several 
AU. The masers are approaching saturation. The magnetic field 
strength is strong (on average 250 mG and as high as 650 mG) 
and shows direction reversal on scales of ~10 12 cm. This can 
be due to the fact that the magnetic field is frozen into a dense 
and turbulent medium although the linear polarization vectors 
indicating the magnetic field direction follow the disk and the 
magnetic field angle with respect to the line-of-sight 9 is ap- 
proximately equal to the disk inclination. The high magnetic 
field strengths indicate that the field is enhanced by a nearby 
magnetic dynamo. 

Field II: This field consists of a newly discovered maser fil- 
ament at -690 AU East of Cepheus A HW2 with a PA= 66°. 
It is likely the result of the shock interaction between the 
HW2 outflow and the surrouding molecular cloud and implies 
a large opening angle (115°) of the outflow. The maser beam- 
ing angle in Field II is ~10" 2 sr with maser path lengths of 
~2 AU while the masers are unsaturated. We find a clear veloc- 
ity and magnetic field orientation gradient along the filament 
consistent with the interaction between a radial magnetic field 



14 



W.H.T. Vlemmings et al.: Magnetic Fields in Cepheus A 



in the HW2 outflow and the magnetic field in the surround- 
ing Cepheus A complex which is almost perpendicular to the 
Galactic plane. The magnetic field strength of 50-70 mG is 
typical for H2O masers found in SFRs. 

Field III: The masers of Field III make up a small part of 
the shell structure found in T01 and, even though our maser ref- 
erence position has an estimated error of up to 25 mas, are fully 
consistent with the shell expansion model parameters estimated 
in G03. We find a magnetic field strength between 30-100 mG, 
consistent with other SFR H2O maser polarization measure- 
ments and find that the magnetic field direction is along the 
shell expansion direction, radial from the central embedded 
proto-star. As these maser have the highest measured linear po- 
larization, Pi = 10%, we can conclude that they are saturated. 
The beaming angle is consistent with a spherical maser geom- 
etry. 

Field IV: Located close to HW2, the maser in this field are 
weak and no magnetic field strength was determined. The up- 
per limits of B\\ ~ 500 mG. The linear polarization indicates 
that the magnetic field is either aligned with the H2O maser 
disk around HW2 or radial toward HW2. 

6. Conclusions 

Strong magnetic fields of up to ~600 mG have been measured 
in the H2O masers around Cepheus A HW2. The strongest 
magnetic fields were measured in the maser structure that was 
identified as a circumstellar disk (G03), suggesting the nearby 
presence of a dynamo source. The field strengths determined in 
the maser regions further from the central source HW2 are 30- 
100 mG, consistent with earlier VLA, VLBA and single dish 
measurements of SFRs. The high magnetic field strengths indi- 
cate that the magnetic pressure is similar to the dynamic pres- 
sure in the outflows around HW2. Thus, the magnetic fields 
likely play a large role in supporting the molecular cloud and 
shaping the outflows in this very active high-mass star-forming 
region. 
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Appendix A: Polarization Modeling and Analysis 

Here we describe the modeling and analysis of the 22 GHz 
H2O maser linear and circular polarization used in this paper 
to determine the magnetic field strength, saturation level and 
intrinsic thermal line width of the maser features. 

A. 1 . Circular Polarization 

For the analysis of the circular polarization spectra we used 
the full radiative transfer non-LTE interpretation, which was 
thoroughly described in V02. There the coupled equations of 
state for the 99 magnetic substate s of the thre e domi nant hyper- 
fine components from Nedoluha & Watsonl ll 19921) (hereafter 
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Fig. A.l. The angle 8 between the maser propagation direction 
and the magnetic field vs. the fractional linear polarization Pi 
for different values of emerging maser brightness (for [r+r v ] = 
1 s" 1 )- The thick solid line denotes the theoretical limit from 
iGoldreich et al.l i 19731) for a fully saturated maser. The shaded 
area is the region of emerging brightness temperatures found 
for the masers in Field II. 
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Fig. A.2. ^-dependence of Eq. IA.1I for increasing emerging 
brightness temperature 7bAQ (for [T + F v ] = 1 s _1 ) from V02. 
The lines for T\,AD. = 10 7 and 10 8 coincide and are the same 
as the lines for lower brightness temperatures. For fully unsat- 
urated masers the dependence is equal to cos 



W.H.T. Vlemmings et al.: Magnetic Fields in Cepheus A 



15 



NW92) were solved for a linear maser in the presence of a 
magnetic field. The emerging maser flux densities of the result- 
ing spectra are expressed in T^AQ., where 7\, is the brightness 
temperature and AQ. is the beaming solid angle. It was found 
in NW92 that the emerging brightness temperature scaled lin- 
early with (r + F v ), which are the maser decay rate F and 
cross-relaxation rate T v . For the 22 GHz H2O masers, F is typ- 
ically assumed to be < 1 s -1 . In star-forming regions it has 
been found t hat T v *2 s~' for T -400 K and F v «5 s 1 for 
T -1000 K JAnderson & Watsoni fl993l) and thus the models 
from V02 (where (r + F v ) =1 s~') have been adjusted to these 
values. 

The model results further depend on the intrinsic ther- 
mal line-width Av th in the maser region, where Av,h - 
0.5(r/100) 1/2 with T the temperature of the masing gas. Model 
spectra for a grid of Av t h between 0.8 and 2.5 km s _1 , corre- 
sponding to temperatures between 250 and 2500 K, were di- 
rectly fitted to the observed I and V spectra using a least square 
fitting routine. As described in V02 the spectral fitting for the 
non-LTE analysis requires the removal of the scaled down total 
power spectrum from the V-spectrum to correct for small resid- 
ual gain errors between the right- and left-polarized antenna 
feeds. This was typically found to be ~ 0.5% of the total power. 
The best fit model thus produced the line of sight magnetic 
field Z?n and the thermal line-width Av t h as well as the maser 
emerging brightness temperature T\,AQ.. However, the uncer- 
tainties in Av t h and T\,AQ. are la rge, as they are str ongly affected 
by m aser velocity gradients dVlemmings & van Langev elde. 
2005f V05). Additionally, T\,AQ. depends on the actual value 
of (F + F v ). We estimate the uncertainties in the fit for Av^ to 
be ~0.3 km s _1 and the uncertainty in log(7bAQ) to be -0.4. As 
the magnetic field depends on the intrinsic thermal line-width 
and emerging brightness temperature (V02), this leads to an 
added uncertainty in the magnetic field determination of -15% 
which has been included in the formal fitting errors. 

When a direct model fit was not possible, we used the re- 
lation between the magnetic field strength B and percentage of 
circular polarization Py. 

P\ — (^max ^min)/Anax 

= 2 ■ Ap-pi ■ fi [G au.s S ]Cos#/Av L [km s -1 ]. (A.l) 

Here 9 is the angle between the maser propagation direction 
and the magnetic field (0° < 9 < 90°) and Av'l is the maser full 
width half maximum (FWHM). V max and V m ; n are the max- 
imum and minimum of the circular polarization and 7 max is 
the maximum total intensity maser flux density. The coeffi- 
cient Ap_p describes the relation between the circular polar- 
ization and the magnetic field strength for a transition between 
a high (F) and low (F') rotational energy level. Ap_p depends 
on Av t h and maser saturation level as described in NW92 and 
V02 as well as on velocity and magnetic field gradients along 
the maser path as shown in V05. We used Ap-p = 0.012, which 
is the typical value we found for the maser of Cepheus A HW2. 
For maser features where no circular polarization was de- 
tected the 3cr upper limits were determined using Eg. IA. II with 
Pv = 6cry// max , with cry being the rms noise on the maser 
V-spectrum determined after Hanning smoothing the spectrum 
(cry -5-8 mjy). 



For maser brightness temperatures 7b AQ > 10 9 K sr it was 
shown in NW92 that the co&9 dependence of Eg. IA. II breaks 
down introducing a more complex dependen ce on 9. This was 
later shown in more detail in Watso n & Wv"klll200ll) for masing 
involving angular momentum J =1-0 and J =2-1 transitions. 
In V02, figure 7 shows the derived magnetic field strength de- 
pendence on 9 toAf_f/ for the 22 GHz J =6-5 transition. This 
figure is repeated here as Fig. IA.2l 

A.2. Linear Polarization 

Maser theory has shown that the percentage of linear polariza- 
tion Pi of H2O masers depends on the degree of saturation and 
the angle 9 between the maser propagatio n direction and the 
magnetic field (e.g. NW92; iDeguchi & Watsoni Il990l) . Figure 
7 of NW92 shows the relationship between 9 and P/ while their 
Figure 8 shows the P; dependence on saturation level. Fig. IA.ll 
shows the dependence of Pi on 9 for various emerging bright- 
ness temperatures as calculated from the models of NW92 and 
V02. A high linear polarization fraction (Pi > 5%) can only be 
produced when the maser is saturated. Additionally, the polar- 
ization vectors are either perpendicular or parallel to the mag- 
netic field lines, depending on 9. When 9 > 6> cr i t ^55° the 
polarization vectors are perpendicula r to the magnetic field , 
and when 9 < 9 cr i, they are parallel jGoldreich et all Il973l) . 
Consequently, the linear polarization vectors can flip 90° at 
very small scal es as was observed in for in stance circumstel- 
lar SiO masers dKemball & Diamond H997I) . 
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